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Abstract 

We study the decays of the J/ if) and ip' mesons to t:~^it~tt^ using data sam- 
ples at both resonances collected with the BES III detector in 2009. We 
measure the corresponding branching fractions with unprecedented preci- 
sion and provide mass spectra and Dalitz plots. The branching fraction for 
J/ijj — > vr+vr^Tr'^ is determined to be 

(2.137 ±0.004 {stat.)tQfSi {sys^-)-QmQ {norm.)) x 10"^^ 

and the branching fraction for ip' — > vr^vr^vr'^ is measured as 

(2.14 ±0.03 {stat.)tQ^^ (syst.)t[!:[]^ [norm.)) x 10~^. 

The J/if} decay is found to be dominated by an intermediate p{770) state, 
whereas the ip' decay is dominated by di-pion masses around 2.2 GeV/c^, 
leading to strikingly different Dalitz distributions. 

Keywords: 

BES III, Hadronic charmonium decays 



1. Introduction 

Previous studies of J/iIj — > tt^tt^tt^ [i|, Q, 0, 0] and V'' — ^ vr+Tr^vr'' 
[l|, 0, ^ found not only an unexpectedly low branching fraction in the case 
of the (world averages: BF{J/^ — > vr+vr-Tr") = (2.07±0.12) x 10"^ and 
BF{il;' — > 7r+7r-7r°) = (1.68±0.26) x 10"^ 0) but also a completely differ- 
ent shape of the di-pion mass spectrum and the Dalitz plot. The fact that 
the /9(770)7r decays as a fraction of all hadronic decays are suppressed by two 
orders of magnitude in the ip' with regards to the J/ if) is especially difficult 
to explain and known as the pn puzzle. Suggested solutions include intrinsic 
charm in the light vector mesons [7,], formation of three-gluon intermediate 



^Some authors [3 claim that the ^' — > tt+tt h'^ branching fraction is as expected and 
the J/i/) — i- 7r^7r~7r'' fraction is much higher than expected for a cc model of the J/i^. 
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resonances [3], a hybrid nature of the ip' [lO|, an additional hadronic ampH- 
tude for the ip' decays the J/ip being dominantly a higher Fock-state 
and so on. 

In this letter we present new measurements of the J/tp — > tt^tt^tt^ and 
ip' — > vr+Tr^vr'^ branching fractions with unprecedented precision using the 
large data samples collected with the BES III detector at the J/ip and ip' 
resonances. These measurements are an important first step to an exper- 
imental inquiry into the puzzle of the decay dynamics, preparing the way 
for a detailed analysis. The large branching fraction of J/ip — > ■K~^'7T~Tr^ 
also makes it an important background process for many other studies (e.g. 
the scalar meson spectrum in J/ip jTr^TT~), an improved knowledge of 
this branching fraction will thus also enhance the precision of those mea- 
surements. 



2. Detector and Monte Carlo Simulation 

BEPC II is a double-ring e+e~ collider designed to provide a peak lu- 
minosity of 10^^ cm'^^s"^ at a beam current of 0.93 A. The BES III [3] 
detector has a geometrical acceptance of 93% of in and has four main com- 
ponents: (1) A small-cell, helium-based (40% He, 60% CsHg) main drift 
chamber (MDC) with 43 layers providing an average single-hit resolution of 
135 //m, charged-particle momentum resolution in a 1 T magnetic field of 
0.5% at 1 GeV/c, and a dE/dx resolution that is better than 6%. (2) An 
electromagnetic calorimeter (EMC) consisting of 6240 CsI(Tl) crystals in a 
cylindrical structure (barrel) and two end caps. The energy resolution at 
1.0 GeV/c is 2.5% (5%) in the barrel (end caps), and the position resolution 
is 6 mm (9 mm) in the barrel (end caps). (3) A time-of-flight system (TOF) 
constructed of 5-cm-thick plastic scintillators, with 176 detectors of 2.4 m 
length in two layers in the barrel and 96 fan-shaped detectors in the end 
caps. The barrel (end cap) time resolution of 80 ps (110 ps) provides a 2a 
K/tt separation for momenta up to ~ 1.0 GeV/c. (4) The muon system 
(MUG) consists of 1000 m^ of Resistive Plate Chambers (RPGs) m nme 
barrel and eight end cap layers and provides 2 cm position resolution. 

For the events to be read-out, one out of seven trigger conditions based 
on combinations of signals from the MDG, TOF and EMG had to be fulfilled. 
At least one of these rather loose conditions should always be fulfilled for 
the events under study, and indeed overall trigger efficiencies very close to 



100% were found for hadronic events containing charged particles 14l |. 



The efficiencies of the detector and the event selection are estimated us- 



ing a Monte Garlo (MG) simulation based on Geant4 [l5|,ll6|. evtgen [17 | 
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is used to generate events; for the J /if) — > 7r'*"7r~7r'^ decay, /3(770)7r events 
give a good description of the data, while for -0' — )■ vr+Tr^vr'^ a mixture 
of p(770)7r and P-wave phase-space events is used d. In both cases, differ- 
ences between the generated and observed distributions are taken care of by 
reweighting the MC events to the data distribution in the Dahtz variables. 
For the estimation of backgrounds, inclusive MC samples are generated by 
KKMC 18,13] — known decays of the J/ip and are modeled by evtgen 
according to the branching fractions provided by the Particle Data Group 
(PDG) [8], and the remaining unknown decay modes are generated with 
LuNDCHARM [17]. Backgrounds from the process J/ij: ^ "yir^TT^ have been 
modeled using amplitudes extracted from a partial wave analysis of BES III 
data. 



3. Data Samples and Event Selection 



This analysis uses a sample of 2.25 x 10^ J/il^ decays [20(] and 1.06 x 10^ 
V'' decays |2lj collected by BES III in 2009. 

Charged particle tracks in BES III are reconstructed using MDC hits. 
We require tracks to pass within ±10 cm from the interaction point in the 
beam direction, within 1 cm of the beam line in the plane perpendicular to 
the beam and to have a polar angle in the range | cos??| < 0.93. Events are 
required to contain exactly one track of positive and one of negative charge, 
corresponding to the 7r+ and 7r~ . Electromagnetic showers are reconstructed 
by clustering EMC crystal energies [22]. The energy deposit in nearby TOF 
counters is included in order to improve the reconstruction efficiency and 
energy resolution. Showers identified as photon candidates must satisfy 
fiducial, timing and shower-quality requirements. Showers from the barrel 
region (|cost?| < 0.8) are required to have an energy above 25 MeV, while 
those in the end caps (0.86 < |cosi?| < 0.92) must have at least 50 MeV. 
Showers from the transition region between barrel and end cap are excluded 
from the analysis, showers within 10° from any charged track. Events 

are required to contain at least two showers fulfilling these criteria. 

For every pair of photon candidates, a full event kinematic fit with the 



^The use of P-wave phase space (and subsequent reweighting in the Dalitz-variables) 
is motivated by the angular distributions and the fact that not much is known about the 
dynamics leading to the accumulation of events around a di-pion mass of 2.2 GeV/c^. This 
procedure has been checked using both toy MC samples and a sample generated using the 
amplitudes extracted from a phenomenological fit to the data. The maximum difference 
in efficiency obtained is taken as a systematic error. 
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x10^ 





Figure 1: Kinematical distributions of the reconstructed vr^s; Top for the J/il) — > ■k'^-k~'k^ 
analysis, bottom for the ?/)' — > 7r+7r~7r° analysis; left showing the tt" momentum, right 
showing the tt" polar angle distributions. 



initial particle (J/V' or ip') four- momentum as a constraint is performed. 
The pair leading to the smallest is kept as the vr'^ candidate if < 50, 
otherwise, the event is rejected. The fit is repeated once more with the 
assumption that the charged particles are kaons; if this leads to a smaller 
X^, the event is also rejected. Yet another kinematic fit is performed with 
the mass of the vr'' as an additional constraint; the resulting x^ is required 
to be less than 50. The invariant mass of the two photon candidates has to 
be compatible with the mass of the 7r°, 0.11 GeV/c^ < m^^ < 0.15 GeV/c^. 

For the V'' — ^ i:~^'k~'k^ analysis, additional requirements are needed to 
suppress backgrounds from radiative decays to e^e~ , /.i^/U" and the J /if) 
and Xc states, namely the invariant mass of the charged pion candidates is 
required to be less than 3 GeV/c^, the energy deposits associated to the 
tracks is required to be less than 0.8 GeV and the penetration depth into 
the muon system less than 40 cm. 
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4. Efficiency correction 

In order to study differences between the simulation and data in track 
reconstruction, an analysis of a specially selected J/V' to Svr candidate sub- 
sample with one or two tracks and two or more photons and with tight 
requirements on one track and the vr'' is performed. Specifically, the TOF 
and ^ information of the charged track are combined to form particle iden- 
tification confidence levels for the vr and K hypotheses; the likelihood for the 
vr hypothesis is then required to be larger than the likelihood for the K hy- 
pothesis. Kinematic fits to the vr'^ mass are performed for all pairs of photon 
candidates, and the pair with the lowest is taken as the vr*^ candidate, if 
the is less than 20. The invariant mass of the object recoiling against the 
system of the track and the TfO is required to be between and 0.2 GeV/c , 
and the recoil direction must be within the tracker acceptance. Using these 
tagged events, the efficiency for finding and correctly reconstructing the sec- 
ond track is determined. The simulation is then corrected as a function 
of polar angle and track momentum to reflect the efficiency found in data 
(which is on average about 2% lower than the simulated efficiency). 

Similarly, an analysis using a subsample similar to the one above but 
requiring two tracks and with tight requirements on the two tracks and 
photons, chosen with the standard photon selection, is performed to test the 
vr" reconstruction efficiency. The tracks are again required to pass particle 
identification requirements as above and in addition are required to have an 
associated energy deposit in the calorimeter of less than 0.8 times the beam 
energy to remove electrons and a penetration depth in the muon system 
less than 40 cm to remove muons. Their opening angle is required to be 
less than 170°. The momentum of the system recoiling against the tracks 
must be larger than 100 MeV/c, and the invariant mass of the tracks and 
the two photon candidates must be above 3.0 GeV/c^. Here the efficiency 
differences found between selected data and simulated events (of the order 
of 0.5%) are used to correct the simulation as a function of vr^ momentum. 
Figured] shows the reconstructed kinematics of the vr*^ for the selected events 
compared to the corrected MC simulation. 



5. Results 

1,859,771 events from the J/ip sample and 7872 events from the ip' sam- 
ple survive all selection criteria. The branching fractions are calculated as 
follows: 

j^p -"set ^ ^ continuum resonance ^ys^ 

■ eMC ■ (-trig ■ -BF(7rO 77) ' 
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Quantity 



Nsel 

^BG 
' resonance 

CMC 



1859771 ± 1364 
8811 ± 1582 
9919 ± 463 
(225.2 ±2.8) Million i 
38.66 ± 0.05 
(100 - 0.2) 



77) 



7872 ± 89 
820 ± 55 
101 ± 32 
106.4 ± 4) Million 
30.91 ± 0.14 
% [14] 



(98.823 ± 0.023)% [8] 



Table 1: Numbers used in the branching fraction calculation 



where Nsei is the number of selected events, N/^^^^^^^^ the number of back- 
ground events from the continuum (estimated from data samples taken at 
3.08 GeV and 3.65 GeV), N^f^^^^^^ the number of background events from 
other resonance processes (estimated using inclusive MC samples) and A^^ 
the number of J/ip or ip' mesons in the sample. eMC is the efficiency de- 
termined from signal MC, etrig is the trigger efficiency (found to be very 
close to 100%, [iJ]), and the branching fraction for vr'^ 77 is taken from 
the PDG — the corresponding numbers can be found in table [TJ In this 
calculation, interference of resonance and continuum processes is neglected. 

Possible systematic errors resulting from the following sources were stud- 
ied: 

• The uncertainty due to the simulation model was estimated by the 
difference in the efficiency with and without the reweighting described 
in section [2] for the J/ip sample and by comparing with the efficiency 
obtained from a sample generated using amplitudes extracted from a 
phenomenological filo for the tp' sample. For both cases, the model 
error is not the dominant systematic error. 

• The absolute energy scale of the EMC is known to an accuracy of 0.4% 

• The photon detection and reconstruction efficiency is described by the 
simulation to within 1% per photon [23]. 



In this fit, contributions from p(770), a higher p with a mass of 2285 MeV/c^ and a 
width of 950 MeV/c^ and a with a mass of 1750 MeV/c^ and a width of 650 MeV/c^ 
were found to lead to an adequate description of the data. 



9 





T /oh 


-)■ 7r~'~7r vr^ 


1p — ) 






Upward 


Downward 


Upward 


Downward 


OUUiCt; Ui bj oLclliciLlC 


^iiaiige 


^iiaiige 


^iiaiige 


v^nange 








(,/oj 


[/(,) 


iViv^ DliliUidLlUli 


u.zo 


n 9*? 
-u.zo 


1 9n 


1 9n 






n 09 

-u.uz 


n 1 S 

U. lO 


n 1 ^ 

-U. lO 




2.04 




2.04 
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Tr*^ kinematic fit 


0.28 


-0.27 


0.27 


-0.27 


tracking efficiency 


1.64 


-1.59 


1.80 


-1.75 


Muon cut 






1.28 


-0.75 


Trigger efficiency 


0.20 


0.00 


0.20 


0.00 


Resonance background 


0.67 


-0.67 


1.45 


-1.45 


Syst. w/o normalization 


2.74 


-2.64 


3.57 


3.33 


Normalization 


1.26 


-1.23 


4.17 


-3.85 


Total syst. uncertainty 


3.01 


-2.91 


5.49 


5.09 


Syst. + stat. uncertainty 


3.02 


-2.91 


5.72 


5.34 



Table 2: Impact of the systematic uncertainties on the measured branching fractions; 
the various sources of systematic uncertainties lead to the listed upward and downward 
changes in the branching fractions. 



• The uncertainty due to the vr finding and kinematic fitting was es- 
timated by performing a different analysis with the J /tp data sample 
(see section |4]) . A tighter selection was applied to the charged tracks 
and no vr'' was required. The difference between data and simulation 
is taken as the systematic error. 

• The uncertainty due to charged particle track finding and kinematic 
fitting was estimated using an analysis with tight requirements on the 
i:^ and one charged track. The efficiencies for finding and reconstruct- 
ing the other track were compared between data and simulation (see 
section H]) . 

• The efficiency of the muon rejection (used only in the ip' analysis) 
was estimated by either dropping the requirement or demanding a 
penetration less than 30 cm instead of less than 40 cm. 

• The trigger efficiency was changed from 100% to 99.8%, reflecting the 
statistical uncertainty of the efficiency determination |14l |. 

• The background from continuum processes was estimated using sam- 
ples taken off-resonance (282 nb~^ of luminosity taken at a center-of- 
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1 2 3 

m(Ti7t) [GeV/c^] 



2 4 6 8 

m(7tV)^ [GeVVd] 




Figure 2: tttt invariant mass distribution (left) and Dalitz plot (right) with backgrounds 
subtracted and corrected for efficiency. Top and bottom graphs show the results for the 
J/xp — > ■K^TT~'K^ and t/'' — ^ 7r"^7r~7r'' analysis, respectively. 



mass energy of 3.08 GeV, compared to 81 pb^^ at the J/^ resonance 
and 43 pb~^ taken at a center-of-mass energy of 3.650 GeV, compared 
to 163 pb~^ at the ip' resonance). The small samples due to the clean 
selection lead to relatively large statistical errors for the continuum 
contribution (18.0% for the J/^ and 6.7% for the if)'). Compared to 
these errors the systematic errors from the luminosity measurements 
or varying beam conditions can be neglected. 

The accuracy of the inclusive simulation for describing background 
from resonant processes was checked in analyses requiring one photon 
less (a TT~^ TT~ 7 final state) or one photon more (a vr"*" vr" vr^ 7 final 
state) and was found to be mediocre; it is assigned an uncertainty of 
100%. 

The normalization (number of J/ij) or tp' events) has an uncertainty 



of 1.23% for the J/V' sample [20t] and 4% for the V' sample [21 1. 



Table [2] shows the impact of the systematic errors on the measured branching 
fractions. 
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The branching fraction for J/ifj — > vr+vr vr'^ is determined to be 

(2.137 ±0.004 {stat.)tlfr^ (sys^Oiaoi {norm.)) x 10"^ 

and the branching fraction for — v vr+vr^vr" is measured as 

(2.14 ±0.03 {stat.)tQ^^ {syst.)^^^^ {norm.)) x 10"^. 

Invariant mass spectra and Dalitz plots are shown in figure [2j The decay 
J/ifj — > -k^tt'tt^ is dominated by p(770) production; the absence of events 
in the center of the Dahtz plot points to negatively interfering higher p states. 
In the case of the -0' — ^ -k^tt^-k^ decay, a small /9(770) contribution can be 
discerned. Most of the events are however clustering around 2.2 GeV/c^ in 
di-pion mass. To disentangle the contributions of various excited p states to 
this peak will require a partial wave analysis. 



6. Conclusion 

The branching fractions for J/0 — > 7r+7r~7r'^ and ifj' — )■ vr+vr^Tr^ have 
been measured with unprecedented precision at the BES III experiment. 
The measurement for ,J /ifj — > t:~^t^~tt^ is in good agreement with the world 
average of BF{J/ip — > tt+tt~7t^) = (2.07 ±0.12) x 10"^ ^] while the result 
for ip' — > TT^TT^TT^ is slightly larger than the world average of BF{ijj' — > 
7r+7r~7r°) = (1.68 ± 0.26) x 10"^ [s^. The ratio of these two branching 
fractions 



5F(^' — > vr+vr-vrO 



BFi j/i. vr+.-vrO) = (^'OO ± O-O^ {syst.))%, 

where correlations between the systematic errors of the two analyses have 
been taken into account, is found to be an order of magnitude smaller than 
the ratio of 12% naively expected from the fraction of decays via three gluon 
exchange. 

The decay dynamics of the J/ip are dominated by the p{770) meson. 
While the p(770) is also visible in the case of the ip' decay, the dynamics 
there is dominated by states at higher masses. Understanding the nature 
of these higher mass states and why they are suppressed in J/ tp decays and 
enhanced in ^' decays may be clarified in a partial wave analysis, which is 
beyond the scope of this letter. 
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